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ABSTRACT: This work investigates the influence of wall treatment on the prediction accuracy of the 
k-e and k-w models, as well as the solid volume fraction, in a solid-liquid fluidized bed. The objective of 
the study is to explore the effects of different y+ values on the performance of these turbulence models 
and the prediction of solid-liquid fluidized bed behavior. Using the Eulerian-Eulerian approach, a 
multiphase flow model was employed to simulate the behavior of solid particles and liquid in a 
fluidized bed. The standard conservation equations of continuity and momentum were solved to 
describe the flow dynamics. The research findings indicate that the k-e model exhibits higher 
sensitivity to near-wall mesh refinement when predicting y+ values compared to the k-a) model. The 
k-w model, however, demonstrates reasonably accurate predictions within the y+ range of 30 to 45. 
Furthermore, the k-e model outperforms the k-w model in predicting solid volume fractions when the 
near-wall region is refined. These results highlight the significance of wall treatment and near-wall 
mesh refinement in accurately predicting turbulent flow behavior in liquid-solid fluidized beds. The 
study contributes to a better understanding of the performance of turbulence models and their 
applicability to solid liquid fluidized beds. 
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1.0 INTRODUCTION 

Turbulence is a common feature of any 
fluidized bed set-up. It causes a chaotic and 
irregular motion of the particles within the 
fluidized medium. It is characterized by the 
formation of eddies, swirling patterns, and 
random fluctuations in velocity and pressure. 
When turbulent flows are bounded by walls, it 
usually presents an additional measurement 
challenge when compared to those in free shear 
turbulence. This results in serious design 
limitations and implementation. Wall effect has 
been found to play a big role on the minimum 
fluidization velocity of fluidized beds. 
Turbulence has also been found to affect 
minimum fluidization velocity [1,2]. There are 
other inherent challenges associated with the 
flow dynamics such as the steepness of the 
mean velocity gradient around the surface and 
the length and time scales of the turbulence 
local to the near-wall [1].The estimation of the 
wall adhesion forces was carried out by 
Wassermnan [3], he compared them with 


hydrodynamic forces. The findings show that 
fluidization characteristics at the boundary of 
micro-flow due to the interplay between the 
ratio of surface and hydrodynamic forces and 
wall effects affect the minimum fluidization 
velocity. Liu [2] studied the minimum velocity 
of fluidization and minimum bubbling velocity 
of silica sand particles in air-blown micro 
fluidized beds in order to appreciate the wall 
effect in micro fluidized beds and, the 
operability and usability of micro fluidized beds. 
The study proposed a better combination of bed 
diameter and particle size in the range of the 
static bed heights from twenty (20) to fifty 
(50) millimeters for the so-called micro 
fluidized bed reactor devised to perform 
reaction(s) with minimal stress from external 
gas mixing and diffusion of gases. 


Shukla et al. [4] investigated the effects of the 
dimensionality of a model (3D or 2D), flow 
regime (laminar or turbulent), grid mesh 
resolutions, and drag law on the void fraction of 
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fluidized beds. 2D simulations were found to 
overestimate the experimental data in terms of 
void fraction distribution across the flow, while 
3D models showed that numerically simulated 
profiles of the void fraction inside the bed 
conform more accurately to experimentally 
derived values. A detailed study of myriad 
numerical model setups proved that for the 2D 
simulations of fluidized beds, there is a mix of 
model parameters such as specularity and 
restitution coefficients, | numerical grid 
resolution, and _ discretization scheme for 
convective terms that tend to produce a good 
correlation between numerical and 
experimental results. Chang et al. [5] found that 
3D simulations exhibited more sensitivity to the 
coefficient of restitution and the specularity 
coefficient than 2D simulations. However, the 
fluidization process in the 2D models was found 
to develop more quickly than that in the 3D 
simulation. 


The effect of varying Reynolds number on 
flow regime of liquid fluidized beds has been 
studied Yao et al [6]. For particle Reynolds 
number less than 20, the porosity is 
comparatively low and the particle dynamics 
are mostly dominated by _ inter-particle 
collisions that bring about anisotropic particle 
velocity variations. When particle Reynolds 
number is increased further to values greater 
than 40, the dynamics of the particles are 
dominated by hydrodynamic effects, which 
leads to decreasing and more anisotropic 
particle variation in velocity. A marked increase 
in the anisotropy happens when the Reynolds 
number of particles increases from 40 to 50, 
and this corresponds to a transition from a 
regime in which collision and hydrodynamic 
effects are equally pronounced to a 
hydrodynamic-dominated regime. With the use 
of the Eulerian-Eulerian two-fluid approach the 
hydrodynamics and heat transfer coefficient on 
various immersed bodies in a fluidized bed 
were studied Fathahi et al [7]. The kinetic 
theory of granular flow (KTGF) was adopted to 
analyze the solid phase. The result showed that 
the spherical immersed body, known to possess 
a better aerodynamic shape, yielded a higher 
heat transfer coefficient in bubbling fluidized 
bed with Geldart B particles classifications. The 
investigation of the hydrodynamic 
characteristics and fluidization behaviour of 
coarse coal particles in a three-dimensional 
liquid-solid fluidized bed operating in near 
turbulent condition has been carried out [8]. 
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The influence of particle size and particle 
density on the homogeneous or heterogeneous 
fluidization behaviour was investigated. The 
findings showed that less dense and fine 
particles are conveniently fluidized, showing a 
certain range of homogeneous expansion 
characteristics, while in the large and heavy 
particles, inhomogeneity may occur throughout 
the bed, including water voids and velocity 
fluctuations. 


Daryus [9] studied the effect of kinetic Prandtl 
number on the k-e turbulence model transport 
equation on the behaviour of fluidized bed flow. 
Factors like the pressure gradient, dissipation 
rate, effective viscosity of the gas and the 
particles in the bed were investigated. The 
investigation observed that turbulent model 
with k-Prandtl of 0.9 gives rise to the most 
accurate results at the superficial velocity range 
of 0.40 m/s - 0.70 m/s, whereas k-Prandtl of 1.1 
produces the most exact results at the 
superficial velocity range of 0.80 m/s - 0.90 
m/s. It was also observed that the decrease in 
the k-Prandtl number leads to a decrease in the 
rate of dissipation; which was the same fact 
with effective viscosity of the gas. There was no 
pronounced difference in the particle’s effective 
viscosity with the change of k-Prandtl number. 
The hydrodynamic behaviour in three-phase 
turbulent fluidized bed contactor with reverse 
current air flow and non-Newtonian liquid was 
observed and juxtaposed with that of 
Newtonian liquid under the same 
environmental conditions [10]. The aqueous 
medium used was carboxy methyl cellulose and 
the hydrodynamic variables studied were bed 
pressure gradient, minimum _ fluidization 
velocity, liquid holdup, bed expansion, and gas 
holdup. It was reported that increasing carboxy 
methyl cellulose concentration increased the 
net drop in pressure through the bed and the 
liquid holdup, while the gas holdup and bed 
expansion were found to decrease. Simulation 
of the heat transfer and flow rate through a 
turbulent gas-solid fluidized bed using the 
Eulerian-Eulerian granular flow model was 
done [11] and it accurately predicted the 
chamber pressure difference, temperature 
change, and bed expansion. 


A comparative study of the flow behaviour in 
three-dimensional supercritical water fluidized 
beds, cold gas-solid fluidized beds, and high- 
temperature gas-solid fluidized beds was 
carried out [12] employing the two-fluid model 
simulations across different fluid velocities. The 
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work also studied particle dynamics, fluid 
velocities, and bed expansions. The outcome 
showed that, the drag forces of particulate 
phase are about 0.4 times particulate phase 
weights in supercritical water fluidized beds. In 
contrast with supercritical water fluidized beds, 
cold gas-solid fluidized beds have lower 
minimum fluidization velocities, the size of 
bubble and slug were smaller, more bubble and 
slug numbers, higher bed expansions and slug 
onset heights, and more uniform §gas-solid 
interactions. There have works on a three- 
dimensional simulation to obtain the voidage 
distributions through a_ time averaged 
simulation. The impact of superficial inlet 
velocity, temperature, pressure, and particle 
size on the radial and axial voidage, probability 
density of voidage, and bed expansion were 
investigated [13]. 

Abdulrahman et al., [14] adopted the 
Eulerian-Eulerian granular multiphase flow 
technique to study the performance of liquid- 
solid fluidized beds. The study quantified the 
characteristics of both the solid and liquid 
holdups. The findings reveal that smaller glass 
beads that represented the solid particles 
require low minimum fluidization velocity; also, 
this velocity requirement increased as the 
diameter of the bead increased. It was also 
found that the bed expansion ratio is directly 
proportional to the velocity of the liquid, and 
inversely proportional to the diameter of the 
beads. Whereas the solid holdup decreased with 
velocity, the liquid holdup increased with 
increasing liquid velocity. 


Resolving near-wall regions properly is very 
important in capturing the effects of the 
different sub layers close to the wall. Wall 
treatment refers to the modeling approach used 
to simulate the behaviour of fluid flow near 
solid boundaries, such as walls, in numerical 
simulations. Accurately representing the flow 
near walls is essential because it influences 
various phenomena, including boundary layer 
development, wall shear stress, and heat 
transfer. In this work wall resolution is done by 
discretization of the computational domain with 
mesh of different sizes and is used alongside 
two different turbulent models, the k—-#£and 


k-w models to study the effect of wall 


treatment and turbulence models on the wall y+ 
values and volume fractions of the bed. 


2.0 METHODOLOGY AND MATERIALS 
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The fluidized bed considered in this work is 
found in [15], as shown in Fig. 1. The continuity, 
momentum and energy equations governing the 
flow in the liquid-solid fluidized bed are solved 
in a 2D bed geometry using ANSYS FLUENT. 
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Fig. 1 Physical set-up of the modeled liquid- 
solid fluidized bed [15] 

From Fig. 1, fluidization occurs in the bed of 
packed particles in a liquid medium. The bed 
section is idealized as a simple 2D geometry. It 
has a column of 0.127m diameter and height of 
1.1m. The liquid is water and the particles are 
glass beads of 1.13mm Sauter mean diameter. 
The model equations were formulated based on 
the following assumptions (i) it is idealized as a 
2D flow (ii) external body forces such as 
magnetic and electrostatic forces are negligible, 
and (iii) the solid particles are of uniform size. 


Using the  Eulerian-Eulerian method 
incorporating the kinetic theory of granular 
flow, the governing equations were formulated 
and shown in Table 1. The closure equations for 
the drag force are presented in Table 2. Two 
types of turbulence models have been studied, 
namely, the standard k-#andk-wtwo- 
equation models for the resolution of the 
transport of kinetic energy and its dissipation 
rate. The standard k — € model assumes that the 
turbulent viscosity is isotropic, that is, the ratio 
of Reynolds stress and mean rate of 
deformation is equal in all directions. For the 
k-—@ model, the first equation takes care of the 


turbulent kinetic energy while the second 
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accounts for the specific rate of turbulent 
energy dissipation. The model constants for the 
two turbulent models are presented in Table 
3,and the fluid properties are shown in Table 4. 
The following boundary conditions were 
applied (i) inlet fluid velocity = 0.076m/s (ii) 
Wall roughness constant is = 0.5 (iii) backflow 
turbulent intensity 


0.001m2/s? (iv) turbulent intensity = 5% (v) 
initial volume fraction = 0 


The equations are solved by the student 
version of the commercial code ANSYS FLUENT 
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2022R. The model equations were solved using 
the turbulent flow option in 2D space and 
applying the second order implicit iterative 
method. An Intel64 bit Operating system 
having Pentium 8 central processing unit 
running on 4GHZ with 4GB RAM was used to 
carry out the simulation. The scalable wall 
treatment was applied while maintaining an 
absolute frame of reference. The iteration was 
run for two different time steps; 10000 and 
20000 at a maximum of 20 iterations per time 
step using atime step size of 0.001secs. 


Table 1: Governing equations 


= (e:p,) +¥.(p,~%))=0 for the liquid phase 


The continuity equations 


= (c,-p,) + ¥.(é,9,v,) =0 for the solid phase 


B (uy — uz) 


The momentum equations 


FE 
2 (e,p,) + V. (e,p,%;) = -5,9P + V.2; + ap.g+ 


2 (z,p, u,)+ V. (z,p,u,) = -2,VP+VP,+V.1, + 


£2, + Blot, - %) 
for the solid phase 


for the liquid phase 


Table 2: Closure equations 
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A, = 3 o2PGosc1 ™ @=)(TF 


Liquid solid drag c,=15 =,(1- =m aoe ExP,|¥, — Vj Gidaspow et al, 1992 [16] 
=a, 
solid pressure P,= a,p,6,+ Zp,(1+ Ps Symlal et al, 1993 [17] 
Radial distribution| | _ |, = 4 Ogawa et al, 1980/18] 
g,=11-C—F] 
may 
solid shear stress | “= Pocet® Pepin * Bey Gidaspow et al, 1992 [16] 
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By ent - 5 P4292)" 2a, 
Symlal1993 [17] 
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Schaeffer, 1987 [19] 
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Bulk viscosity Lun et al, 1984 [20] 
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Table 3: Closure constants for the k —= and k —-# models 


oy = 1.0 
6,=1.30 a= 0.56 
C,=1.44 f=0,3 
C€z,=1.92 gp =0.9 
Table 4: Fluid and solid Properties 
Density 999.5Kg/m3 
Specific heat capacity 1006J/kg.K 
Viscosity 0.001236Kg/m.s 
Molecular weight 28.966Kg/Kmol 
Thermal conductivity 0.0242W/m.K 
Density 2540Kg/m3 
Sauter diameter 1.13e-3 m 
Packed bed volume fraction 0.6 


3. RESULTS AND DISCUSSION 


The results obtained from the study are 
presented and discussed in Figs. 2-8. The y* 
values represent the distance from a 
computational grid point to the nearest wall, 
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normalized by a characteristic length scale in 
the flow. It provides information about the flow 
physics near the wall and is used to determine 
the appropriate grid resolution in the near-wall 
region. 
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Fig. 2 Variation of y+ values with bed height for k-e turbulence model at 2.5mm grid size 
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Fig. 2 shows the variation of yt values against 
bed height for 2.5mm grid size. It can be seen 
that the values increase with bed height and 
reached well over 40 within a short distance of 
the bed height, and thereafter started 
decreasing until it reached 10 at about 0.5m bed 
height. Above 0.5m, the y+ values maintained a 
constant value. This means that the first 
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computational cell centroid is located relatively 
far away from the wall. The y* value of 40 falls 
within the "log-law" region, which is relatively 
far from the wall and outside the viscous 
sublayer. This suggests that a global mesh size 
of 2.5mm has a coarse resolution near the wall 
which might compromise near-wall prediction. 


7.000 8.000 


Wall Yplus (liquid) 


Fig. 3 Variation of yt values with bed height for k-e turbulence model for 2.5mm grid size 
with a refined mesh at the wall 


Fig. 3 shows the effect of near-wall mesh 
refinement on the turbulence. From the wall, 
the y* values increased with bed height until at 
about 0.3m when it became unstable with 
almost unnoticeable marginal increase. The 
maximum y* value was as low as 9.3, indicating 
that near-wall flow behaviour, and resolution of 
the effect of viscosity and velocity gradient can 
be done with high fidelity. 


(mi) 
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For the k-e model with a mesh size of 1.25mm 
shown in Fig. 4, there is no marked uniformity 
in the pattern of variation of the y+ values near 
the wall. However, a maximum value of about 
25 can be observed. This suggests that the first 
computational cell centroid is located within the 
viscous sublayer, and it shows a clear difference 
of over 50% compared to the 2.5mm mesh. 


Fig. 4: Variation of yt values with bed height for k-e turbulence model at 1.25mm grid size 
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The k-w model gave a lower y* value of about 
35 as can be seen in Fig. 5 for the 2.5mm mesh 
size which suggests better resolution of near- 
wall region. When the global mesh size was 
changed to 1.25mm, the y* values for the k-w 
model moved closer to the 30 mark as shown in 
Fig. 6. It is evident that the k-e model showed 
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more sensitivity to near wall mesh refinement 
than the k-w model in predicting the y* values. 
Past authors have suggested different wall y+ 
resolutions for near-wall flows. Reports of 
available work have been compiled and 
compared with the present work as shown in 
Table 5. 
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Fig. 5: Variation of y+ values with bed height for k-w turbulence model at 2.5mm grid size 
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Fig. 6: Variation of yt values with bed height for k-w turbulence model at 1.25 grid mesh 


Table 5: Comparison of different near wall treatments and y* values with present work 


| Muhammad [21] 2019 k-« Standard wall function 6.7 
Salim et al [22] 2009 k-¢ Standard wall function 30< y* < 60 
Saverio et al [23] 2018 not specified Not specified 30< y* < 100 
Present work 2023 k-e« Standard wall function 9.3<y*<45 
k-aw Standard wall function 30< y* < 40 


magnitude of fluctuations of liquid volume 
fraction is higher for both the k-e and k-w 


The vector plot of volume fractions in the bed 
is shown in Figs 7 and 8. It shows that the 
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models near the wall. This is consistent with the 
findings of Marefatallah et al [24] which showed 
that the highest magnitude of fluctuation is 
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usually around the near-wall region and can be 


up to 3 to 4 times greater than that at the center 
of the bed. 


alm 


Fig. 7: Vector plot of volume fraction for k-e model and 2.5mm mesh 
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Fig. 8: Vector plot of volume fraction for k-w model and 2.5mm mesh 


4.0 CONCLUSION 


A numerical simulation of the near-wall 
treatment in a liquid-solid fluidized bed has 
been undertaken. From the two comparative 
simulations carried out involving k-£ and 
k-@ models, it was observed that the 
k-—£model showed more sensitivity to near 
wall mesh refinement than the k- @ model in 
predicting the y* values. The k-w model gives a 
reasonably accurate prediction in the range 
of30< y* < 45. The best prediction of liquid 
volume fraction was made by the k- model 
when the near-wall region was refined. This 
means that different turbulence models have 
different capabilities and promises in the 
prediction of the accuracy of numerical 
simulations of multiphase flows, and the choice 
of a more appropriate model solely depends on 
the particular parameter of interest. This study, 
therefore, establishes the fact the kK-—€ and 
k-—w models are good closure models for two- 


phase flows in liquid-solid fluidized beds. 
However, it is suggested that more turbulence 
models be studied and _ perhaps further 
refinements near the wall be done to establish 
the best model for the fidelity of CFD 
simulations of fluidization process in fluidized 
beds. 


Definition of symbols 


Symbol Meaning 

T Phase stress-strain tensor, 

B Interphase momentum 
exchange coefficient 

Loss Radial distribution function 

ess Coefficient of restitution for particle 
collisions 

0; Granular temperature 

Vr Kinematic viscosity 

€ Volume fraction 

p Density 

V Velocity 

g Gravitational acceleration 

tT Shear viscosity 

a 
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P Pressure 

y Bulk viscosity 

- Dynamic viscosity 
Re Reynolds number 
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